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— We present details of our recent calculation of a 2 corrections to the para- 

i ' positronium decay into two photons. These corrections are rather small and 

j^. our final result for the parapositronium lifetime agrees well with the most 

c~| , recent measurement. Implications for orthopositronium decays are briefly 

discussed. 
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I. INTRODUCTION 

Precision measurements with nonrelativistic bound states provide a sensitive test of 
bound state theory based on Quantum Electrodynamics (QED). Theoretical predictions 
generally agree well with experimental results. However, there are a few observables where 
the agreement is not satisfactory. The best known example is the lifetime of the orthopositro- 
nium (o-Ps), where some experimental results |l|,§] differ from the theory by several standard 
deviations. The theoretical prediction is not yet complete, because two-loop QED correc- 
tions to the o-Ps decay have not been evaluated yet. However, those corrections would have 
to be very large to reconcile the theory with experiment. 

Recently, we have reported a calculation of corrections to the parapositronium (p-Ps) 
decay into two photons @ in the second order in the fine structure constant a. We argued 
there that the smallness of the 0(a 2 ) corrections to the p-Ps lifetime makes it unlikely 
that the analogous effects in o-Ps could alone explain the discrepancy between theory and 
experiment. 

The purpose of the present paper is to present a detailed description of the calculation 
reported in || . We begin with a brief discussion of the current theoretical and experimental 
status of positronium (Ps) decays. 

Theoretical predictions for p-Ps and o-Ps decay rates into 2 and 3 photons, respectively, 
can be expressed as series in a: 
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are the lowest order decay widths of the p-Ps and o-Ps, respectively, and the ellipses denote 
unknown higher order terms which we will neglect in our analysis. Corrections of 0(a) 
were calculated in |4[] for p-Ps. For o-Ps the most accurate result was obtained in , where 
references to earlier works can be found. The logarithmic two-loop correction was found 
in for o-Ps and in for p-Ps. The leading logarithmic correction at three loops was 
computed in ||. Some partial results on the 0(a 2 ) corrections for both p-Ps and o-Ps can 
be found in P,|9|-|I"2"f . A complete calculation of B p has been reported recently in but the 



value of B has not been obtained so far. 

Eqs. ([j],|2|,|3|) gi ve the following predictions for the lifetimes: 



a^ 2 



Tf™Y = 7989.42 /xs" 1 + T^B P (-) = (7989.42 + 0.043 B p ) /is~\ (4) 

r^ r / = 7.0382 /is" 1 + T^B (-) = (7.0382 + 0.39 ■ 10~ 4 B ) /is" 1 . (5) 

On the experimental side, the Ann Arbor group MM has found 



T^p^gas measurement) = 7.0514(14) /xs 1 , 
T^pg (vacuum measurement) = 7.0482(16) /is -1 , (6) 

which, for B = 0, differ from (|5]) by 9.4<r and 6.3a respectively. This apparent disagreement 
of experiment with theory has been known as the "orthopositronium lifetime puzzle" for a 
long time. More recently, an independent measurement by the Tokyo group found 



n* p Ps (Si0 2 measurement) = 7.0398(29) fis~\ (7) 

which agrees with the theory if B is not too large.Q Thus, the present experimental situation 
is rather unclear. If future experimental efforts confirm the Ann Arbor results (§), the 
orthopositronium lifetime puzzle could be solved if B turns out to be unusually large, 
e.g. ~ 250 for the vacuum measurement. Alternatively, one might speculate that some 
"New Physics" effects such as o-Ps decays involving axions, millicharged particles, etc., 
cause the excess of the measured decay rate over the QED predictions. Some of those exotic 
scenarios seem to have already been excluded by dedicated experimental studies. (For a 
review and references to original papers see e.g. |i"5| .) 

p-Ps can decay into 2 photons and its lifetime is too short to be measured directly. For 
a long time its precise value remained unknown. However, it was realized in |16[ that such a 



measurement would be very useful, since the calculation of the coefficient B v is much easier 
than that of B . Therefore, p-Ps offers an easier precision test of the bound-state QED. This 



observation motivated the measurement of the p-Ps lifetime [[16] which is 6.5 times more 
accurate than the best previous results. 

To fully utilize that experimental result and enable the rigorous test of bound-state QED 
envisioned in fllEf , we undertook a complete calculation of the 0(a 2 ) corrections to p-Ps rate 



H . We found that the non- logarithmic part of those corrections is small; the coefficient B p 

is 

B p = 5.1(3), (8) 

and the theoretical prediction for the p-Ps lifetime becomes (half of the logarithmic a 3 term 
in (|l|) is taken as an error estimate) 

r t p h !p r s y = 7989.64(2) /is" 1 . (9) 



Comparing this number with the most recent experimental result |]T6 

r^L p P9 = 7990.9(1.7) fis- 1 (10) 

we find excellent agreement between theory and experiment. 

In the next Section we explain our approach to Ps decays with the example of the leading 
order calculation. In Sections [TQ| and |V| we discuss, respectively, the so-called soft and hard 
O (a 2 ) contributions. The final result is presented in Section M. 



1 Another measurement of the o-Ps lifetime was performed earlier in Mainz |TJ|] with the result 
Pq^Ps = 7.031(7) ixs" 1 . This result is consistent with the theoretical predictions but uncertainties 
are clearly larger. For a summary of older experimental studies see |L5|]. 



II. FRAMEWORK OF THE CALCULATION 

Positronium is a bound state of an electron and a positron. Its energy levels and lifetimes 
can be well understood within the framework of nonrelativistic expansion in QED. The decay 
width of p-Ps into two photons can be written as: 
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where P = 
momenta, q 
the process e + e~ — > 27, A denotes polarizations of the photons, and 
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(M p _p s , 0) is the four-momentum of the positronium, k\ and fc 2 are photons' 
(p) is the positronium nonrelativistic wave function, A(X,p) is the amplitude of 
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is the spinor part of the p-Ps wave function. The wave function of the p-Ps is normalized 
to unity, 
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hence the usual normalization factor, l/2M p _p s , is absent in Eq. (|TT|). From Eq. (pT|) one 
sees that r p -p s is determined by two quantities: the annihilation amplitude of the electron- 
positron pair into a pair of photons, and the positronium bound state wave function. Since 
the typical velocity of electron and positron in Ps is O (a), the annihilation amplitude can be 
expanded both in a and in the relative momentum p of the electron and positron. Correc- 
tions to the positronium wave function are computed using the standard time-independent 
perturbation theory and Breit correction to the nonrelativistic Hamiltonian. 

We begin with the calculation of the decay rate to leading order. Since the electron and 
positron in Ps are nonrelativistic, we employ non-covariant perturbation theory. Although 
it is not the only way, the non-covariant technique permits an easier evaluation of the soft 
corrections to the p-Ps decay rate. 

The on-shell amplitude of the process e + e~ — > 77 reads 
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Here E p = y/m 2 + p 2 is the energy of the electron (or positron), A±(p) are the projectors 
on the mass shell, 



A±(P)= 2 [ 1± — K- 



(15) 



and k is the three-momentum of the photon in the final state. 



To the leading order one can neglect the dependence of the annihilation amplitude on 
the momentum |p| ~ ma, which is small in comparison with m and \k\ ~ m. One obtains 
the following leading order amplitude: 

Aiiot 
^lo = -(ae 2 )(afe)(aei). (16) 



From Eqs. (|lq , |12D it follows that the photon polarization vectors are perpendicular to each 
other, i.e. e.\e 2 = 0. 

Computing the trace in Eq. (O) and integrating over p one obtains 



J-^Tt[A lo (\,p)^p}<P(p) 



2^,2 



which gives the following leading order decay width: 
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m 2 
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The calculation of higher order corrections to positronium lifetime is performed in the 
framework of the Nonrelativistic Quantum Electrodynamics (NRQED) || which we regu- 
larize dimensionally ]T7|] . In this framework, all contributions are divided into soft and hard 
corrections. The soft contributions come from the momenta region of the order of k ~ ma 
and thus are sensitive to the details of the bound state dynamics. On the other hand, the 
hard corrections arise as contributions of the relativistic momenta k ~ m; their effect can 
be described by adding <5(r)-like terms to the nonrelativistic Hamiltonian. 

Special care in the present calculation is required because of a single 75 matrix in the 
positronium wave function. Since a consistent treatment of single 75 is known to be a 
problem in dimensional regularization, below we describe how we have dealt with it. 

We use the following fact: if the p-Ps decays into two photons with polarizations e,\ and 
e2, then the photons are polarized so that e,\e2 = 0. Let the first photon be polarized along 
the x and the second photon along the y axis; the three-momentum of the first photon k 
being along the z-axis. If we introduce two auxiliary vectors 

m m 

then the standard four- dimensional representation of the 75 matrix, 75 = ijolxlylz, can be 
written as: 

75 = l^UM- (19) 

We use this equation to define 75 in d dimensions. 

With an additional trick it is possible to avoid any reference to photon polarization 
vectors. Let us consider Eq. ( |IT1) and use there the explicit representation for the 75 matrix, 
Eq. (p!9|). Then the following trace has to be computed: 

Tk[Mii4M*p\ e^el (20) 



Taking into account that e±e2 = e\^k = e\^P = 0, one concludes that the final result can 
only depend on e\ and e\. In this situation, one can average Eq. (|20D over directions of e\ 
and e2, provided one respects the above constraints. We arrive at the following formula: 

TrlMtMMa/,] e a A - ^WlHM^\ Tg, (21) 

where 

(3 - 2e)dgrfg - d^O* - dffi 
"" 4(l-2e)(2-e)(l-e) ' 

^ = 5V ~ ^ ~ %Vv ( 22 ) 

According to the discussion of different momentum regions contributing to the second 
order corrections, we divide the coefficient B p in Eq. (Q) into three parts: 

B p = B s p ^ Icd + 5j ard + ^ oft , (23) 

where i^^ uared is the contribution of the one-loop amplitude squared and £?p ard ' soft are the 
hard and soft two- loop contributions. The square of the one-loop amplitude is easily obtained 
from the one-loop result: 




^squared / = j gggg ^ 



In contrast to the first order correction which arises from hard photon exchange only, 
the second order correction is more difficult to compute, because of the appearance of the 
soft scale effects. In addition, hard and soft corrections are not finite separately. Below we 
discuss the calculation of these corrections. 



III. SOFT SCALE CONTRIBUTIONS 

As follows from Eq. flTTf ), in order to obtain Bp° ft one has to compute relativistic correc- 
tions to the annihilation amplitude e + e~ — > 77 and relativistic corrections to the positron- 
ium wave function induced by the Breit Hamiltonian. Accordingly, the soft contribution is 
separated into two pieces: 

Bf' = B s p oit (A A) + Bf\WF) . (25) 



A. Relativistic corrections to the amplitude 

The calculation of relativistic corrections to the amplitude is straightforward. One starts 
with the on-shell amplitude Eq. (0) and expands it up to relative order 0{p 2 /m 2 ). The 



calculation of these corrections can be performed in three dimensions. To demonstrate this, 
let us write the correction to the amplitude in the form:0 



,(2) 



SA = AXjViVy (26) 

To calculate the correction to the p-Ps lifetime induced by SA, we have to compute the 
following integral: 

To this end we use the Schrodinger equation in the momentum space: 



d d p , . 9 f d d p Anamp 2 f d d k d>(k) . , 
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Rewriting 
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shifting the integration momenta in the first term p — > p + k and using the fact that the 
scaleless integrals vanish in dimensional regularization, we arrive at 



dp 

(27r) d 



p)p 2 = mEil>(p). (29) 
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We see that the amplitude A {j is needed only in the e — > limit, where it can be easily 



calculated. We obtain: 



3m 2 



SA=-J^A h0 , (30) 



which induces the following 0(ma 2 ) correction to the p-Ps lifetime: 

Bf\AA) = y. (31) 

Recently there has been some discussion in the literature [|ll],|ll],|il|] concerning the lin- 
early divergent integral in Eq. (p9|). Our approach to the linear divergence is based on di- 



mensional regularization which permits a consistent treatment of hard and soft corrections 
simultaneously. We dealt similarly with linearly divergent integrals in our recent calcula- 



tions of the 0(ma ) corrections to positronium S-w&ve energy spectrum (T9|,|2(| and found 
agreement with earlier results obtained in a different regularization scheme [ 21f| . This gives 
us confidence in the result given in Eq. fl2"9"|). 



2 A linear term in the expansion of the amplitude in p does not contribute to the decay rate since 
the positronium ground state wave function is spherically symmetric. 



B. Relativistic corrections to the wave function 

Relativistic corrections to the positronium wave function can be computed using Breit 
Hamiltonian. Since we regularize all divergences dimensionally, we need that Hamiltonian 
in d dimensions, as it has been derived in ||20|| . 

Breit Hamiltonian projected on the S-states reads 

tti \ P 4 d — 1 \p 2 .1 dna .. . ira ir . .... 

U(r ' p) = -^ + ■4S-{m' C(r »} + ^ (r) -4*rf[ ffi '^W'^ (r »- (32) 

Here the Pauli matrices <Tj and er^ act on the two-component spinors of the nonrelativistic 
electron and positron respectively, and 



aV(d/2 

n d/2-l r u 



C(r) - ~^m^ < 33 > 



is the (^-dimensional generalization of the Coulomb potential. The operator U(r, p) induces 
the following correction to the p-Ps width, 

£- = Ai + A a , (34) 

where the first term is due to the 5(r)-part of the operator U and the second term is 
due to the remaining terms of that operator. We note also, that Tlo in the above formula 
stands for the leading order decay width computed in ^-dimensions, in contrast to the three- 
dimensional result Eq. (pf). All necessary formulas can be extracted from the calculation 



described after Eq. (36) of Ref. [|(J. One obtains: 

a 2 n 



Ai = -— (--41n(ma)-2) 



A2 = T~(7~ 41n(mc0 + f)' (35) 

Thus the wave function correction contribution to the positronium lifetime becomes 

Bf\WF) + 2n 2 In - = t. + 2tt 2 In — + ^, (36) 

y a 2e ma 8 

On the LHS of the above equation we have separated the logarithm of the fine structure 
constant to be consistent with division of corrections introduced in Eq. ([I]) . 

C. Final result for the soft contributions 

The sum of the corrections to the annihilation amplitude (|31|) and to the wave function 
(J36l) gives the final result for the soft contributions, 

ir 2 1 07ir 2 

B ^ = --2-K 2 \nm+-^. (37) 

p 2e 24 v ; 



IV. HARD SCALE CONTRIBUTION 

The second class of corrections are the hard scale contributions. The corresponding Feyn- 
man diagrams are shown in Figs. |2|(VP), |3|, |], |5|. They should be computed in dimensional 
regularization, with external electron and positron at rest. 

Bp ard consists of three types of contributions: vacuum polarization insertions in the 
photon propagators, light-by-light scattering diagrams, and two-photon corrections to the 
annihilation amplitude, 

5j ard = 5j ard (VP) + 5j ard (LL) + 5j ard ( 77 ). (38) 

Vacuum polarization insertions into the one-loop graphs (an example is shown in 
Fig. |(VP)) were computed in HHH, 



5£ ard (VP) = 0.4473430(6). (39) 



A. Light-by-light scattering contributions 

One class of the second order corrections to the p-Ps lifetime arises from the photon- 
photon scattering, shown in Fig. |[ These contributions are relatively small; it is interesting, 
however, that the positronium lifetime measurement of increased precision may become 
sensitive to effects of non-linear QED. 

For the diagrams shown in Fig. ^, we find that the planar ones are equal (a=b), and 
so are the non-planar ones (c=d). Further, both classes remain unchanged when we cross 
the internal photons. Therefore we only need to compute two diagrams, one planar and one 
non-planar. We will call them L\ and L 2 . The contribution of those diagrams to the final 
expression for the amplitude is 4Li + 2L 2 . The factor 4 arises because we have two types 
(a) and (b) which differ by the orientation of the fermion loop, and both can have parallel 
or crossed internal photon lines. For the type L 2 we only have a factor 2 for the orientation 
of the fermion loop, because crossing of the internal photon lines simply exchanges (c) and 
(d). 

The light-by-light diagrams in Fig. |3] have an imaginary part because of the two-photon 
cut. This complicates a numerical integration over Feynman parameters and we used a 
different method to evaluate them. The idea is to (formally) assign a large mass M to the 
internal fermion line. The diagram may then be expanded in ratio m/M using the so-called 
Large Mass Expansion P^-F2B| . This reduces the task to the calculation of two-loop vacuum 



or massless propagator integrals. The price to be paid is that the result is a series in m/M, 
while we are interested in the value of the series at m/M = 1. 

Fortunately, using symbolic manipulation programs one can compute many terms of the 
series; in our calculation about twenty terms were computed for each diagram. The resulting 
series converge well, especially for L 2 , and the number of computed terms is sufficient to 
obtain an accurate estimate of this contribution at x — 1. The behavior of the series for L\ 



is improved if we make a change of variables m/M = Jz/(2 — z). The n— th term of the 
series in variable z decreases roughly like In n/n 2 . Finally, we find 

9 



In = 0.494(30), L 2 = -0.348(20), (40) 

and the contribution of these diagrams to the coefficient B p is: 

5j ard (LL) = 4Li + 2L 2 = 1.28(13). (41) 

B. Two-photon corrections 

Another class of corrections is generated by the two-photon diagrams, shown in Fig. |j. 
These are the most difficult diagrams we have to compute, since in general they diverge 
and a regularization is required. There arise two types of divergences: First, there are the 
ultraviolet (UV) divergences; they are removed by an appropriate renormalization. Second, 
there are infrared (threshold) singularities. They will remain in the final result for the hard 
scale contributions and vanish only in the sum with the soft contributions. 

The evaluation of the hard corrections is made possible by a combination of the analytical 
and numerical methods. The idea is to construct an infrared safe expression from divergent 
Feynman amplitudes by subtracting appropriate counterterms, which can be computed an- 
alytically. The construction of the counterterms is based on the following observation: in a 
given Feynman diagram the infrared singularities appear when the loop momenta are small. 
In this situation the propagator of the virtual electron in the t-channel can be contracted to 
a point. As the result, the infrared behavior of those Feynman diagrams is identical to that 
of the two-loop three point functions considered earlier in p7| , |20| . 

We proceed in the following way: after constructing an infrared finite expression, we 
combine propagators using Feynman parameters; perform momentum integrations analyt- 
ically; extract ultraviolet (UV) divergences and integrate numerically over typically 5 (in 
some cases 6) Feynman parameters in the finite expressions. For the numerical integration 
we use the adaptive Monte Carlo routine VEGAS ||28|| . 

Let us illustrate the basic steps of the calculation by considering as an example the non- 
planar box diagram D\ shown in Fig. [|. A power counting shows that this diagram is UV 
finite but IR divergent. To demonstrate how the IR counterterm is constructed we consider 
a symbolic expression for this diagram (after taking the trace over Dirac matrices): 

J (2n) D (2n) D P i mi + 2ph)(q-2pl 2 )(Q + 2ph)(q-2hp)(ll + 2p 1 l 3 + 2m>)' l > 

Here p = (m, 0) is the four momentum of the incoming electron or positron, I3 — h + 1% is 
the sum of the loop momenta l\p, and p± = p — q, where q is the four- momentum of the 
outgoing photon. The quantities f^' in the numerator denote the uniform functions of the 
loop momenta: 

/ (0 (AJi,Afe) = A7 (<) (Ji,J2). (43) 

Only terms with f^ and f^ diverge in IR. We use the following identity: 

D 1 = Di i ^ + (Dt°' 1) -\Dt 0,1) ] ) + M W) l • (44) 
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where the counterterm 



D\ 



(i=0,l) 



ct 



is obtained by expanding the propagator of the electron 



in the t-channel in Taylor series in small loop momenta, 



D\ 



(<=o,i) 



d D h d D l 2 



f(0) + /( D 



ct 2m 2 J (2tt) d (2tt) d l\l\(l\ + 2 V l 1 ){\\ - 2pl 2 )(lj + 2pl 3 )(l 2 " 2Z 3 p) 

x (l - £±^| . (45) 



2m 2 



Examining Eq. ( ]4"4| ) one recognizes that the first two terms in that equation are finite, both 
in the UV and the IR, and hence can be evaluated numerically. The last term, D± ' , 

L J ct 

is divergent. Since the t-channel propagator has been contracted to a point, this term 
corresponds to a three-point, rather than four-point Feynman amplitude. Such integrals 
were computed in a previous study 



the counterterm 



D\ 



(j=o,i) 



see also fl20| ). Using those results one can obtain 



ct. 



analytically: 



D 



(i=0,l) 



ct £ 



4 + 48 + 12vr 2 In 2 + 6tt 2 + 42C 3 . 



(46) 



Similar procedure was applied to evaluate the remaining Feynman diagrams Di. In some 
cases, like e.g. for the planar box diagram D 2 , the overall subtraction is not sufficient and a 
more sophisticated approach is required. The results of the calculation are summarized in 
Table 0. 

Finally, we have to consider the one-loop diagrams. We need the results including terms 
(9(e), because they will be multiplied by divergent renormalization constants. The results 
are summarized in Table |I[ 

For the renormalization one needs the electron wave function renormalization constant, 
SZ e = Z e — 1, and the mass counterterm 8m computed in dimensional regularization to 
0(a 2 ). These results can be found in [29]. For completeness, we collect here the relevant 
formulas: 



8Z e 

5m 

m 



a 6Z® 



a 2 8Z® 



m 



+ a z 



Sm^ 



m 



(47) 



where 



szP 



SmW 



m 



(4tt)^/ : 

-2-4 

e 



-771 



(-26) 



T + 
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Smf® 



9 45 _ 

2? + 47" 



49tt 2 

4 
17tt 2 



+ 16vr 2 ln2 



„ . * 433 

24C 3 + --, 



+ 8tt 2 In 2 - 12C 3 + 



199 



m 2e* 4e 4 

The final result is obtained by putting the many pieces together: 



(48) 
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m 4 ^ d ( TT ) = 5Z^B + — B 1+ ( S -^—\ 2 B 2 
F m \ m J 

+8ZV> (s 1 + S 2 + S 3 + — B ± ) 
\ ml 



J') 



"' \i=5,7,ll,12,15,17,19 / i=l 



TV 2 



— -42.19(27). (49) 

Here 5j are the tree-level diagrams (Fig. P, Si are one-loop diagrams (Table |), C$ are the 
one-loop diagrams with mass insertions (Table [TI|), and D^ are the two- loop diagrams (Table 

0)- 

For the complete hard correction we add Eqs. ([^, ^l|, ^9|) and find 



TV 2 
? hard _ " , o_2 



B™ Td = + 27r 2 lnm- 40.46(30). (50) 



V. FINAL RESULT 

The final result for the second order non-logarithmic correction to the p-Ps decay rate 
into two photons is obtained as a sum of the soft ( |37|) and hard fl50D pieces, and the square 



of one- loop corrections ( 24]) . Adding them one finds 

B p = 5.1(3), (51) 

and the theoretical prediction for the p-Ps lifetimes becomes 

r* h _!£f = 7989.64(2) ^s" 1 . (52) 

In this equation we have not included the contribution of the decay mode p-Ps— > 47. This 
decay channel increases r p _p s by approximately 0.01 /is -1 p0|j31[| . 

VI. CONCLUSION 

We have described the numerical and analytical methods employed in our calculation of 
the second order QED corrections to the p-Ps lifetime. Our final result, the non-logarithmic 
correction 0(a 2 /tt 2 ) with the coefficient B p = 5.1(3) increases the decay rate by approxi- 
mately 0.1 yus -1 . The resulting theoretical prediction is in excellent agreement with experi- 
ment. 

The framework of this calculation is the Nonrelativistic QED with dimensional regular- 
ization. The dimensional regularization facilitates the separation of scales, the cornerstone 
of the effective theory. Its main technical advantage is that no additional scales, such as a 
photon mass, are introduced by the regularization, so that only single scale integrals need 

12 



to be computed. Unfortunately, these integrals are still too complicated to be evaluated 
analytically. We computed them numerically by subtracting IR counterterms. Those were 
constructed using simpler integrals which are known analytically, as we explained in Sec- 



tion [TV]. The multidimensional finite integrals were computed using adaptive Monte Carlo 



integration routine VEGAS [28 



The approach described in this paper can be also applied to the calculation of the 0(a 2 ) 
corrections to the o-Ps decay into three photons. In particular, the IR counterterms can be 
constructed in a similar manner. However, the larger number of diagrams and two additional 
integrations over the three-photon phase space make this problem significantly more difficult. 

Somewhat surprising is the smallness of the second order corrections found in this paper. 
It resulted from a strong cancellation between soft and hard pieces computed in dimensional 
regularization. We would like to stress that the soft and hard pieces are not separately finite 
and depend on the regularization (in this sense they are "scheme-dependent"). For this 
reason, large constants accompanying divergent pieces in Eqs. ([5l],[J7]) may have no direct 
physical meaning. Unambiguous information is provided by the scheme-independent results, 
Eqs. (p^j31|j39|J^l| ) , which are all of order [several units] x (a/n) 2 . 



In the absence of a complete result on 0{a ) to r o _p s it is interesting to discuss what 
our result for p-Ps might imply for the o-Ps lifetime puzzle. Although nothing can be said 
rigorously, we believe that our result indicates that no dramatic enhancement in the 0(a 2 ) 
effects in o-Ps decay is possible. In Table fV] we have summarized available results on the 
second order effects for both p-Ps and o-Ps decays (we have not included the partial results 
for the soft corrections in o-Ps since they are scheme-dependent). One can see that, with 
the exception of S squar , all radiative corrections are comparable for o-Ps and p-Ps decays. 
On the other hand, significantly larger value of B squa,rcd for o-Ps can be traced back to 
a larger value of the one-loop correction to o-Ps— > 37 rate. The relation of the one-loop 
corrections for p-Ps and o-Ps decay, however, is rather natural since the number of diagrams 
is approximately three times larger for o-Ps decay. Thus, apart from the factor related to 
the number of Feynman diagrams, there seems to be no significant difference in the structure 
of radiative corrections to o-Ps and p-Ps decays. 

Therefore, it is difficult to imagine that a complete calculation of the 0(a 2 ) correction to 
o-Ps decay will result in a dramatically large number, necessary to resolve the o-Ps lifetime 
puzzle. We believe that this puzzle will be solved by continuing experimental studies and 
we look forward to learning their results. 
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FIG. 1. Tree-level amplitude of the p-Ps decay, higher order mass counterterms, and their 
values. 
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FIG. 2. One- loop corrections to p-Ps decay. 
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FIG. 3. Light-by-light scattering contributions to p-Ps decay. 
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FIG. 4. Two-loop photonic diagrams for the p-Ps decay. 
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FIG. 5. Mass counterterms in one-loop diagrams. The numbering conforms to Fi£ 
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TABLE I. Values of one-loop diagrams 
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TABLE II. Values of counterterm diagrams 
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TABLE III. Values of two-photon diagrams in Fig. f|. 
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TABLE IV. Comparison of available results for p-Ps and o-Ps. For the various coefficients B 
we use the notation introduced in the text. B n ' 1 describes multiphoton processes, p-Ps — ► 47 and 
o-Ps — ► 57, and B d (LL) is due to effects of non-linear QED: 7*7* — ► 77 for p-Ps and 7* — ► 37 
for o-Ps. Tlo denotes the lowest order decay width of a given state. 
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